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Abstract 

A frequency-domain procedure for the design of sliding mode controllers for multi-input, 
multi-output (MIMO) systems is presented. The methodoK)gy accommodates the effects of 
parasitic dynamics such as those introduced by unmodeled actuators through the introduction of 
multiple asymptotic observers and model reference hedging. The design procedure includes a 
frequency domain approach to specify the sliding manifold, the observer eigenvalues, and the 
hedge model. The procedure is applied to the development of a flight control system for a linear 
model of the Innovative Control Effector (ICE) fighter aircraft. The stability and performance 
robustness of the resulting design is demonstrated through the introduction of significant 
degradation in the control effector actuators and variation in vehicle dynamics. 

Introduction 

A reconfigurable control system is one which is able to compensate for sudden, 
potentially large, unknown failure o^ents in real-time using on-line adaptive control laws and/or 
adaptive redistribution of control effort. The objective of reconfiguration is system stability 
white retaining some level of required performance and handling qualities. As a research area, 
reconfigurable flight control has seen rapid growth in the past decade. The motivation for 
developing reconfigurable flight controls is clear. Failures during flight are inevitable, especially 
in combat aircraft. If the flight control system is capable of stabilizing the aircraft and providing 
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acceptable handling qualities (if the aircraft is piloted), it may be possible to save the airframe 
and save lives. Many reconfigurable algorithms that have been discussed in the literature 
involve, to some degree, (1) failure detection, (2) system identification, and (3) flight control law 
reconfiguration.''^ The time necessary to complete these operations can be critical, especially if 
the airframe in question is open-loop unstable. Even direct adaptive methods like dynamic 
inversion with neural nets to remove inversion error*''^ require some convergence time for 
adaptation. The method of sliding mode control (SMC) offers a characteristic which makes it 
very attractive for a reconfigurable application. The variable structure nature of the SMC 
controller allows it to adapt to very large parameter variations instantaneously. In fact, an SMC 
controller is invariant to so-called matched uncertainty. Plant uncertainties are defined as 
matched when they lie in the image of the plant input matrix, i.e., the uncertainties affect the 
plant dynamics only through the plant input channels. Thus, if the system is invariant in the 
presence of uncertainties such as those arising from airframe damage, there exists no need to 
perform failure detection, system identification, or on-line control law redesign. Indeed, it was 
this potential simplicity and instantaneous adaptation time that motivated the research to be 
described. 

The most significant shortcoming with SMC designs in control applications is their 
inability to accommodate unmodeled parasitic dynamics. In flight control applications, these 
dynamics are typically those of the actuators driving the control effectors. One method which is 
known to help alleviate this problem is the use of asymptotic observers.'^ A previous work by 
the authors demonstrates an observer-based SMC for a SISO flight control application and 
introduces a frequency domain approach for designing the sliding manifold. This approach has 
been extended to a MIMO observer-based SMC controller for an aircraft with stable 
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unaugmented dynamics.'^ In the current work, the unstable highly-coupled vehicle dynamics of 
the tailless Innovative Control Effector (ICE) fighter aircraft are addressed, and a frequency 
domain procedure for the important selection of observer eigenvalues is presented. In addition, a 
method of observer loop shaping which is equivalent to model reference hedging is introduced. 

MIMO Sliding Mode Control 

Overview of SMC 

There are several excellent survey articles regarding SMC theory and applications. 

Hence, only a brief, simplified overview is given here, with emphasis on implementation and 

design issues. The key properties of a sliding mode controller are well known and are reviewed 
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here without proof 

a. ) While on the sliding mode, the system dynamics are invariant to matched uncertainty. 

b. ) The hypersurface that describes a = 0 defines the transient response of the system 

during the sliding mode. 

c. ) While on the sliding mode, the trajectory dynamics are of a lower order than the 

original model. 

Consider the uncertain system with m inputs and n states given by: 

i(0 = A(x,?) + B(jc, t)M(0+ (1) 

where / represents the parameter uncertainties present in the system and is assumed to be 
unknown but bounded by some known function of time, system state, and control vectors. In 
simple terms, the objective of SMC is to define 

a. ) m sliding surfaces or manifolds, represented in vector for as cr(x) = 0, and 

b. ) a variable structure control given by 

M(x,t) = /7Sgn (o) (2) 
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such that the system is driven to the sliding surface (t=0 in finite time and remains upon this 
surface for all subsequent time. When <7-0, a sliding mode is said to have been obtained. 
Equation (2) is said to describe variable structure control because the control structure is 
dependent upon the sign of the function cr. This problem statement implies a two-step design 
process. First, the sliding manifold(s) must be designed. This can be accomplished by a wide 
variety of approaches ranging from arbitrarily selecting desired error dynamics to an LQR-like 
design approach utilizing the state equations in the so-called regular form. A method which 
combines a square system feedback linearization approach ’ ' with traditional frequency 
domain loop shaping is used in this work. The second step in the design process is to select p 
such that the reaching condition is met. This can be done analytically using a Lyapunov stability 
criteria or simply by use of computational simulation. 

SMC Implementation Issues 

The control law in Equation (2) is undefined while <t- 0; therefore, the control effort 
required to maintain the ideal sliding mode is discontinuous with an infinite frequency switching. 
This control signal is unacceptable for actual mechanical systems and a continuous 
implementation of a SMC is almost always sought. One of the most common approaches to 
obtain a continuous control signal is the inclusion of a so-called boundary layer near the sliding 
surface. This can done by replacing the signum function in Equation (2) with a finite-slope line 
with limits at ±1 as shown in Equation (3). 


m(x, t) = /7 sat 


^ a 

K E 


( 3 ) 


Another implementation issue with SMC controllers is the problem of unmodeled 
parasitic dynamics. The interaction of these parasitic dynamics with the SMC controller (even 
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when a boundary layer is used) results in chatter and very often causes system instability. A 
simple solution to the actuator problem might be the inclusion of the actuator in the SMC design. 
Unfortunately, this is difficult to implement in practice because of the increased order of the 
sliding manifold associated with the additional dynamics. In general, the order of the manifold 
will increase by the same order as the modeled actuator dynamics. This means, for a second 
order actuator, at least two derivatives of the output signal are required. For a real system with 
measurement noise, these additional derivatives make this approach unattractive. There are 
several other approaches which have been proposed.'^ One which is both effective and 
intuitively appealing is the use of an asymptotic observer. 

Observer-Based SMC 

Asymptotic observers construct a “high-frequency bypass loop,” essentially hiding the 
unmodeled parasitic dynamics from the SMC controller. Selection of the observer eigenvalues 
is a crucial part of the observer-based SMC design. If the observer eigenvalues (and hence 
gains) are too large, chatter and instability results. In addition, sensor noise becomes a problem. 
If the observer eigenvalues (and gains) are too small, robustness to system parameter variation is 
lost. A method is sought to select the observer gains. Since a MIMO square system architecture 
effectively decouples the control v^;jables, it can be shown that multiple observers (one for each 
control variable, and each with its own eigenvalues) are desired for reducing the interaction of 
each sliding mode control action with the parasitic dynamics. Such an architecture is shown in 
Fig. 1 and is directed toward the specific flight control application to be discussed in the design 
example. Consider the following system definitions corresponding to Fig. 1 
.ice9t", 
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Linear Plant 

X = A ;c + B 
y = Cy X + Dy^a 

^1 ~ ^zl ^ ®zl ya 
^2 = Cz2 ^ + Dz2 ya 
Zy - Cj3 ®z3 ya 

Reference Model 

X, = A,x, +B^r^ 
y, = C, X, + D, r„ 

Observer 1 

xj = (a qi — G I C^i )x/ + Bqj K + G I 7 / 

V/ = C„i X/ 

Observer 2 

^2 = (Ao 2 “ ^ 2 Cz2 )^2 + Bo 2 K by,. + G 2 2^ 

>^2 = Co 2 ^2 
Observer 3 

^3 = (^o3 ~ ^ 3 Cx 3 )xj + B „3 K + G 3 zj 
h = Co3 ^3 

Output Feedback 

y = yi +h+h 


Actuators 

= AaXfl + Ba KfeJ'c 
ya ~ ^a ^ ®a^b 


Col e with zeros in the rows 

corresponding to states not output by 
Observer 1 . Gi = Observer gains 

Co 2 e with zeros in the rows 

corresponding to states not output by 
Observer 2. 62 = Observer gains 

Co 3 e with zeros in the rows 

corresponding to states not output by 
Observer 3. G 3 = Observer gains 


Figure 2 shows an equivalent systei^ with unity feedback. It is relatively simple to derive the 
state space representation for the system in Fig. 2 as 


Xe = AeXe + Beyc 


y= 




-[“ySl 


yc 


y=CcXe 


X,. = 


X 

Xa 

^1 


(5) 


where 
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It is instructive to now examine the transfer functions :^-(s) and compare them with both the 

>c 


original nominal plant transfer functions with no actuator and the plant with the actuator. Figure 

3 shows the Bode plots of -^(s)for the ICE vehicle, where y a = angle of attack. The nominal 

Yc 

plant transfer function is shown along with the observed feedback signal for various observer 
bandwidths. Several observations are worth noting. First ol all, while the addition of actuator 
dynamics (with no observer) changes the apparent relative order for the SMC, the apparent 
relative order of the observed signal- is the same as the nominal plant with no actuators. This is 
important because the sliding manifold is one degree less than the relative order of the plant and 
is not tolerant to changes in relative order. Second, fast observers (those with eigenvalues far 

into the left half plane) display a characteristic distortion in the Bode magnitude curve of ^(s) 


near the frequency of the actuator dynamics and a large increase in phase lag. Third, by slowing 
down the observer, the distortion in the magnitude curve and the large phase lag can be greatly 
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reduced. This fact can be used to help determine the best observer gains. Similar Bode plots are 
examined for each controlled variable in the square system architecture. 

Model Reference Hedging 

Observer-based SMC helps deal with the problem of unmodeled parasitic dynamics. 
However, the addition of the observer does not directly address the issue of control saturation. 
The observer does help with rate saturation to some degree because of phase lag reductions at 
medium to high frequencies. However, position saturation is still an issue especially when the 
system does not have redundant control effectors. Initially, in an attempt to deal with control 
saturation, a method called model reference hedging was pursued. The concept of hedging has 
been successfully demonstrated in a dynamic inversion design approach.^^ The concept in words 
is this: “The reference model is moved backwards (hedged) by an estimate of the amount the 
plant did not move due to system characteristics the control designer does not want the adaptive 
control element to ‘know’ about. The actual accelerations are subtracted from the expected 
accelerations (assuming no actuators). This difference represents the amount of desired 
acceleration which was not achieved due to the actuators and should capture nonlinear saturation 
of the actual actuators. This difference is then subtracted from the reference model acceleration. 
Since all this is done in a dynamic^nversion setting, these accelerations are pseudo-commands 

for the dynamic inversion controller, and Johnson et al. call this “Pseudo-Control Hedging.” The 

26 

concept has great merit and is very successful in such applications. 

In order to employ this method in an SMC design, some modification to the 
implementation is required. Rather than subtracting the acceleration difference from the 
reference model pseudo-command, the controlled variable output is subtracted from the expected 
output (through a nominal system with no actuators). This signal (yh) represents the amount of 
unachieved performance due to the presence of the actuator. The signal yh passes through a 
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hedge gain (Kh) and is subtracted directly from the second order model reference states as shown 
in Fig. 4. 

In order to investigate how hedging affects the system, an analytical expression for the 
system with hedging is needed. Consider, initially, a model architecture as shown in Fig. 5. In 
this system, the hedge plant, Gh, is simply the nominal plant model (assuming it is stable). The 
subsystem block labeled “Hedged Reference Model” is the subsystem shown in Fig. 4.. Since 
hedging occurs on individual control variable channels, it is possible to write the transfer 
function for the hedged reference model. It can be shown that this transfer function (for a 2"^* 
order reference model) is given by 


Xrh = 


mo 


2 

s +m|S + m 


oj 


2s +m|S 


s +mis + m 


o; 


Kh Yh 


( 6 ) 


= Gr(s) Xc -Gf(s) Kh yh 


Equation (6) indicates that the hedged reference signal consists of two parts, a reference model 
part and a hedge signal passing through a “hedge filter,” G|. With this definition, the block 
diagram can be redrawn as shown in Fig. 6. With the following hedge system definitions 


Hedge Plant Hedge Filter 

+B, 

Ty, X, +D, yj =C,x^ + H,K^y^ 


Output Feedback 

h = h + h + h + yf 


( 7 ) 


the equivalent unity feedback system (as in Fig. 2) is defined in state space as 
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( 8 ) 


BD, Kb 
Ba Kb 

B„i Kb + G, D^i D, Kb 
Be= B„2 Kb + G2D,2 Ky, 

B„jKb + GjD,3D, Kb 

Bf Kb Db Kb 

BbKb 

Ce=Oniyxn Bniyxn^ C^j C^2 ^o3 

De = DfKhDhK^ 

Now the effects of hedging can be examined on the Bode plots of the transfer functions 
— (s) as was done with the observer. After examining several different systems, it was noted 
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that the hedge signal has basically the same loop shape lor each system. It resembles a 
derivative at low frequencies; it peaks, and then rolls off at the relative order of the plant at high 
frequencies. This is not surprising considering the hedge filter has the same form in all cases — it 
resembles a high-pass filter. All real plants have a high frequency roll-off. When these two are 
placed in series, the result is a characteristic “hump” shape, fhe question then is this; can the 
same beneficial effects of hedging be achieved with a simplified “equivalent hedge” transfer 
function? This is highly desirable because, in its current form, hedging can only be used with a 
stable plant. 

Consider a hedge filter of the form (a high-pass filter): 

Gf(s) = ^- (9) 

s + af 

And a hedge plant of the form (a low-pass filter): 

Gh(s)= 3 — (10) 

s +ahs + bh 

The hedge plant, Gh, in Eqn. (10) has relative order 2 and is intended to be used with a system 
with a relative order of 1 . The extra pole is added in order to have additional high frequency 
roll-off In general, the simplified^model can be designed to be close to the original hedge 
system and will achieve the same desired effects. Note, the actuator states are not required as 
they are in pseudo-command hedging.^^ The following technique is proposed for creating the 

yh 

hedge model. Begin by plotting the Bode plots of — (*) and the nominal system as before with 

yc 

zero hedge gain. The basic form of the desired hedge model loop shape is as follows: 

(1 ) -1-20 dB/dec slope at low frequencies 
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(2) -20 r dB/dec slope at frequencies where the actuators distort the magnitude curve 
(r = relative degree of nominal system with no actuators) 

(3) -20 r-20 dB/dec slope at high frequencies 

The desired hedge signal loop shape can be created (for a system with relative degree 1) using 
the hedge plant and hedge filter given in Eqns. (9) and (10). The pole in the hedge filter is 
placed at the high frequency end of the magnitude distortion and the two poles of the hedge plant 
are placed at the low frequency end of the distortion. A conceptual example of this is illustrated 
in Fig. 7. With the correct loop shape, the hedge gain is then adjusted to minimize error in the 
phase plot and the distortion in the magnitude plot. If the observer is very fast, or if no observer 
is used, there is a minimum hedge gain that will stabilize the system. Using this method of 

examining the Bode plots of — (^) , it is easy to determine the minimum hedge gain. Examining 

yc 

the Bode plots will also give an indication of the upper limit for this gain. In general, if the 

Yh 

hedge signal moves the magnitude plot of (^) above the nominal plant magnitude plot, the 

yc 

closed-loop system will have unacceptable overshoot and phase lag. This can lead to instability 

if outer control loops are closed around the SMC system. 

A 

In MIMO applications, the cross-coupled transfer functions need to be examined as well. 
If cross-coupling effects are strong, the interaction of the command for one channel with the 
unmodeled parasitic dynamics can lead to instability in another channel. Therefore, if necessary, 
a hedge model for the cross-coupling term can be designed and the feedback loop properly 
shaped. The same guidelines introduced above are also used for the cross-coupled hedge 
models. This cross-term hedge signal is then added to the primary hedge signal. 
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It is worth noting that the proposed hedging method no longer resembles the concept of 
hedging given in the initial reference work.^^ Note in Fig. 1 1 that all the hedge signal dynamics 
occur in parallel with the observer loop and no longer enter the reference model. The hedging 
method here is really a form of observer loop shaping. In fac t, this approach is similar to Loop 
Transfer Recovery (LTR) as used in Linear Quadratic Gaussian (LQG) control in which 
optimality of the observer is traded for increased stability margins. The hedge signal used here 
also has the effect of “tuning down” the observer at certain frequencies in order to recover 
desirable stability margins. It can also be shown that the hedge system is attempting to invert the 
actuator dynamics. This is why it is highly dependent on the actuator bandwidths. When viewed 
in this light, one could argue that this is simply another form of an SMC prefi Iter-type design.’^ 
The preceding discussion has attempted to illustrate the effectiveness of model reference 
hedging and has demonstrated a design procedure— all in the frequency domain. Observers and 
hedging in the feedback paths allow the SMC to see the plant (nearly) without parasitic 
dynamics. Thus, some of the performance and robustness of the SMC design approach are 
recovered. Unfortunately, rigorous proofs that quantify stability and robustness for nonlinear 
applications are, as yet, unavailable. Some initial observations about robustness using a 
frequency domain approach are fortfiScoming.^^ 

Frequency Domain SMC Design Procedure 

The following technique is based upon the design approach offered in previous 

applications.'"*'^ 

(1) The vehicle model is obtained, along with an estimate of the frequency beyond which 
parasitic dynamics (or unstructured uncertainties) are likely to come into play. This 
frequency is referred to as the limit frequency in this discussion. Actuator dynamics are not 
included in the nominal plant model. 
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(2) A reference model is chosen for each control variable channel. Since piloted flight control is 
of interest in the present application, this reference model should be selected with an eye 
towards Level I handling qualities with no pilot-induced oscillation (PIO) tendencies. This 
can be accomplished based upon a pilot model-based handling qualities and PIO prediction 
technique.^*’^^ 

(3) The desired feedback structure of the control system is determined with a square system 
architecture. For example, if a roll-rate command flight control system is desired, then roll 
rate (pc) becomes the output of the reference model, and estimated roll rate ( p ) is fed back to 

the SMC system from the observer. System error is then defined as e(t) = P(,(t)-p(t). 

Special care must be given to ensure the dynamics of the uncontrolled variables remain 
stable. This assumption is analogous to the minimum phase assumption used in feedback 
linearization.^^ Also, if there are redundant control effectors for the desired 
moments/pseudo-commands, a control distribution matrix must be defined. 

(4) The sliding manifold, a, is chosen based upon the following principles: 

a) a is derived from a tracking error expression as 

a = e (t) V + K ^_,e (t)"-^ + K„e (t) + K , J e (t) dt (11) 

where p is the relative order of the system, i.e., the number of times the vehicle output 
must be differentiated for the input to appear. Note that the (p-1)^* derivative of the error 
signal is included in the definition of ct. An integral term also appears in Eqn. (11) to 
counter the steady-state bias often created with the use of a boundary layer. 
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b) Recognizing that a boundary layer is to be implemented, the control law is 
expressed as a linear transfer function, assuming the boundary layer is large enough to 
remain within the linear region of the saturation element. 


u(s)=-o 

s 


K ^ 


^K^sP-'+K sP-^+K + 


p- 2 


e(s) 


( 12 ) 


The parameters Kj are chosen to provide desirable properties in the frequency domain. 
This means creating a loop transmission with broad K/s-like characteristics around 
crossover.^*^ This will always be possible since enough derivatives are included in 
Eqn. (12) to create exact K/s characteristics beyond a certain frequency (at least as high 
as the limit frequency). Parasitic dynamics are deliberately excluded in this formulation. 
This step will involve obtaining an estimate of Kp, as this value will determine the 
crossover frequency of the loop transmission. This crossover frequency is selected to 
provide acceptable stability margins as obtained from a Bode plot of the loop 
transmission, but using a value of Kp at least as large as the largest amplitude limit of any 
of the control effectors. Thd latter criterion is included to accommodate maximum trim 
positions of the control effectors. As opposed to typical designs involving loop shaping, 
very high crossover frequencies may result from this step. Indeed these frequencies may 
be well beyond the limit frequency. This result is of no immediate concern. If a MIMO 
system is being designed, a classic sequential loop closure technique is used, thereby 
sequentially and independently determining the coefficients for each sliding manifold. 
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(5) Using the Kj’s just determined in the definition of the sliding function, the existence of a 
sliding mode is verified in the inner loops using a true SMC. This step is completed without 
the observer, actuators, reference model or pilot model, i.e., assuming that no outer-loop is 
being utilized. If necessary, p is increased until sliding behavior is created. The initial value 
of p = Kp obtained in step (4b) should be considered a lower limit in this process. While an 
analytical approach to determine p is certainly possible here, a more expedient route of 
establishing the sliding mode using a computer simulation of the system is also possible. 
Near perfect tracking in the face of large parameter variations should be observed. The 
control signal, however, will exhibit very high frequency switching. 

(6) A boundary layer is included in the controller by replacing the signum function sign(a) with 

P 

While maintaining an approximate constant - = Kp, 

E 

increase the boundary layer thickness, e, until no high frequency switching is evident. 
Again, a simulation of the SMC system is a convenient way of finding this e. Near-perfect 
tracking (with a continuous control signal) in the face of large parameter variations should 
be observed. 

A 

(7) Parasitic dynamics are included in the model. The SMC controller will very likely be 
unstable at this juncture. 

(8) An asymptotic observer is created for each control variable channel as discussed above. 

(9) A hedge model is designed as described above. 


the saturation function sat 
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(10) The frequency domain characteristics of the open and closed-loop SMC system with 
observer, boundary layer and reference model are examined to ensure that stability of the 
linear system is in evidence. 


A criticism that can be levied against SMC design procedures such as the one discussed here 
is that, after introduction of a non-adaptive boundary layer, one has essentially created little more 
than a MIMO PID controller. However, an examination of the characteristics of the equivalent 
unity feedback systems proves this assumption incorrect. Cimsider again the SMC system of 

y 

Fig. 1 . Focusing upon any one of the closed-loop transfer func tions of the MIMO system, — (s) , 

r 


one can create an equivalent unity feedback loop transmission as 

T _ (y/*‘) 
i-(y/r) 


(13) 


This can then be used to create an equivalent serial compensation element, (s), which 
represents the combination of the controller and the observ^er. 



(14) 


By maintaining the SMC and observer parameters constant, but introducing changes in the 
vehicle or actuator dynamics, it can be shown that (s) will vary in a manner that attempts to 
compensate for changes in the vehicle dynamics. This behavior is attributable to the presence of 
the observer(s) in the feedback loops. Of course, the changes in G^^ (s) are instantaneous. 

Design Example: ICE MIMO Linear 6-DOF Model 

ICE Vehicle 

The Integrated Control Effectors (ICE) aircraft model has been developed by Lockheed 
under an Air Force Research Laboratories (AFRL) sponsored program^' and is the vehicle of 
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choice for many controls applications in the current literature.'^ ' It is a single engine, 
multi-role, supersonic, tailless fighter aircraft with a 65 degree sweep delta wing (see Fig. 8). As 
illustrated in Fig. 9, the conventional control effectors include elevens, symmetric pitch flaps, 
and outboard leading edge flaps. The innovative control effectors include pitch and yaw thrust 
vectoring, all moving tips, and spoiler slot deflectors. The all moving tips and spoiler slot 
deflectors have zero lower deflection limits. 

The static aerodynamic force and moment data were collected by NASA Langley 
Research Center and AFRL using wind tunnel tests with a 1/1 S"* scale model. Additional wind 
tunnel tests during Phase II of the ICE program provide uptlated data for simulation models.' 
There are strong multi-axis effects and highly nonlinear interactions between the close-coupled 
control surfaces. The full nonlinear simulation of the ICE vehicle is proprietary; however, 
linearized models at various flight conditions (which do not include nonlinear effector 
interactions) have been made available to this research directly from the NASA Langley 
Research Center. 

SMC Design 

The SMC design considers simultaneous control of longitudinal and lateral degrees of 
freedom. The response variables afe angle of attack, a, roll-rate about the velocity vector, ps, 
and sideslip, p. Decoupled tracking performance is desired. 

(1) The flight condition used for this study is: I g, wings-level, Mach No. = 0.3, 

Alt = 15,000 ft. The linearized model for this flight condition contains eight states, eleven 
outputs, and eleven control inputs. The open-loop plant is unstable and highly coupled. 

(2) Reference models for each response variable arc created as well-damped second- 
order systems with dynamics predicted to yield Level 1 handling qualities with no PIO 
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tendencies. The methodology utilizing a pilot structural model is employed.^* The reference 


models so obtained are: 


G„(s) 

Gp(s) 

Gp(s) 


— (s) = 


a. 


— (s) = 

Pc 
Pr 


100 

(s^ + 20s + 100) 

10 0 

(s^ + 20s +100) 


.. 10 0 

(s' + 25s +100) 


(15) 


(3) A square system architecture is defined as shown in Fig. 10. As shown, the feedback 

structure of the flight control system consists of angle of attack command and hold, roll-rate 

command, and sideslip command and hold. Analytical models of the human pilot controlling the 

angle of attack (a) and roll-attitude (<t>) loops are developed and included in the Simulink 

simulation of the pilot/vehicle system.^* A common crossover frequency of 1 .5 rad/s are chosen 

for the a and (|)-loops. The P loop is not closed by the pilot, i.e., it is assumed that the pilot flies 

the vehicle with “feet on the floor.” The controller to be designed generates pseudo-commands 

for angle of attack, stability axis roll rate, and sideslip angle. These demands are allocated to the 

eleven control effectors using a pseudo-inverse approach.^^ If the preferred actuator positions 

are assumed to be zero, the solution to the control allocation problem using a pseudo-inverse is 

A 

B, =W;'B^(BW„'B'’)-' (16) 


where the weighting matrix, Wu, is assumed to be identity. The resulting control distribution 
matrix, B|„ is an 11x3 matrix with constant elements. Noncorrelated sum-of-sines inputs serve 
as tracking commands to the pilot in the (() and a loops, p command is zero. 

(4) The system (without actuators) has a relative order of 1 for all three control variables. 
Therefore, the form of the sliding manifolds are 
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a = K„e + K^j jedx; e = {(«,-«), (p, - rt (pr " 


(17) 


The control laws, assuming the use of a boundary layer, then can be expressed in linear form, 
resulting in one zero to be placed and a gain to set for each loop during the loop-shaping design. 
Using a traditional sequential loop shaping technique, the parameters of Eqn. (17) are 
determined— with roll rate being the first loop to be closed, followed by a; then p. The designed 
manifolds are given in Eqn. (18), and the design Bode plots for the last of the loop closures 
(P-loop) are shown in Fig. 1 1 . A crossover frequency of co = 1000 rad/s is set for each channel. 


Uc«(s) = 5000 

Ucp(s) = 1000 

Uc/j(s) = 5000 


's + 10 ^ 

, s j 

s + 20 j 
^s-rlO^ 


(18) 


(5) Sliding behavior is verified by simulation. As expected, the SMC provides near- 
perfect tracking and is invariant to the system parameter changes. The control output shows the 
classic high frequency switching. The controller achieves decoupled tracking of a, ps and p, as 
desired. 

A 

(6) Next, the boundary layer is increased until a continuous control signal is achieved. 
For this model, e = 1 is chosen for all three channels. Again, the performance is excellent, even 
in the face of large system failures. 

(7) The next step is the inclusion of the actuators. The actuators are second order with 
rate and position limits. As expected, the nominal system immediately goes unstable. 

(8) -(9) The next steps are to design the observers and hedge models. In this procedure, 
the eigenvalues of the observers are chosen to be real, and nearly identical. Using the frequency 
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domain approach introduced earlier, it is seen that (without hedging) the a-channel observer 
must have eigenvalues of 10 rad/s or less; the p-channel observer can have eigenvalues as large 
as 100 rad/s; and the P-channel must have eigenvalues of 5 rad/s or less. It appears that hedging 
will help in the a-channel, is unnecessary in the p-channel. and is probably required in the 
P-channel. Three cross-coupling hedge models were found to be necessary in the design. The 
parameters chosen for the observer and hedging models are shown in Table 1 . The design Bode 
plots for each primary channel are given in Figs. 12 - 14. 

Simulation Results 

To illustrate the robustness of the observer-based SMC design, the system is exercised 
with measurement noise, system failure, and the pilot model in the loop. The Simulink 
simulation is run with an ODE2 solver using a fixed time step of At = 0.0005 sec. Vehicle 
failure occurs at t = 20 sec. Vehicle failure is defined as: 

(1) Plant failure (A-matrix is multiplied by 2 — with the exception of the elements 
describing kinematic relationships; B-matrix is multiplied by 0.75) 

(2) All actuators experience a 50 ms time delay 

(3) Left elevon rate limits are reduced from 150 deg/s to 10 deg/s; position limits are 

A 

reduced from ±30 deg to ±15 deg 

(4) Symmetric pitch flap undamped natural frequency is reduced from 63 rad/s to 
1 0 rad/s 

(5) Left leading edge flap is jammed at +5 deg 

(6) Pitch nozzle actuator undamped natural frequency is reduced from 39 rad/s to 
1 0 rad/s 
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(7) Yaw nozzle actuator undamped natural frequency is reduced from 39 rad/s to 
1 0 rad/s 

The resulting outer loop tracking is shown in Fig. 15. After system failure, tracking is noticeably 
degraded, but the vehicle remains stable. The actuator deflections are not shown, however, it is 
seen that nearly all the actuators are in nearly constant rate saturation after the failure. 
Predictions of handling qualities and PIO susceptibility for the SMC system after the failure can 
be made using the pilot/vehicle analysis technique.^*'*® As seen in Figs. 16 and 17, Level 1 
handling qualities are still predicted for the a-tracking task, although PIO susceptibility is 
increased. Level 11 handling qualities are predicted for the <|)-tracking task, and PIO tendencies 
are moderately strong. In viewing these results, certain observations are in order. First of all, 
this is a significant system failure, and it is remarkable that the controller is able to maintain 
stability. Second, the input maneuvers are very aggressive. After a failure of this magnitude, a 
pilot would not be commanding such a demanding profile. Third, the HQ and PIO level 
predictions assume a static pilot model. In actuality, a real pilot would compensate for the 
perceived change in vehicle dynamics. Therefore, the results shown above are conservative. 
One can compare these results with a classic loop-shaped controller. As Fig. 18 indicates, this 
design goes unstable shortly after th^ failure occurs. 

Conclusions 

Using asymptotic observers with SMC has been shown to help mitigate the adverse 
effects of parasitic dynamics, and a design method for choosing observer gains has been 
presented. A form of model reference hedging has been shown to be equivalent to an SMC pre- 
filter which helps shape the feedback loop— thus partially removing the effects of actuator 
dynamics. A complete design procedure has been presented which incorporates a frequency 
domain approach to select 1) the sliding manifold, 2) the observer eigenvalues, and 3) the 
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hedging dynamics. The approach was successfully applied to the design of a flight control 
system for a model of the Innovative Control Effector fighter aircraft. Stability and performance 
robustness was demonstrated in the presence of significant actuator degradation and variation in 

vehicle dynamics. 
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Figure Captions 

Figure 1 Multiple observer-based SMC system 
Figure 2 Unity feedback equivalent of Fig. 1 

ct 

Figure 3 Bode plots, — , ICE, various observer eigenvalues 

Uca 

Figure 4 2"^* order hedged reference model 

Figure 5 Hedged system architecture 

Figure 6 Equivalent hedged system architecture 

y 

Figure 7 Creating the hedge model, — 

Yc 

Figure 8 Innovative Control Effector (ICE) vehicle 
Figure 9 ICE control effectors 

Figure 10 Simulink diagram of the ICE pilot/vehicle/SMC system 

Figure 1 1 Loop transmission for p-loop of step (4) in ICE SMC design, all other loops closed 
Figure 12 Bode plots, — ^ , ICE hedge design 

Figure 13 Bode plots, — , ICE hedge design 

^cp 

Figure 14 Bode plots, , ICE hedge design 

Ucp Y 

Figure 15 Outer loop tracking (a, (]), and P), ICE, with SMC. actuators, observers, hedging and 
noise, failure at t = 20 sec 

Figure 1 6 a-tracking task HQ and PIO predictions, ICE failed system 
Figure 17 <|)-tracking task HQ and PIO predictions, ICE failed system 

Figure 1 8 Outer loop tracking (a, <(), and p), ICE, with loop-shaped controller, actuators and 
noise, failure at t = 20 sec 
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Table 1 : Observer and Hedge Models for ICE 
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Figure^; Hedged system architecture 
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Figure 6; Equivalent hedged system architecture 
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Figure 8: Innovative Control Effector (ICE) vehicle 
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Figure 9: ICE control effectors 
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Figure 10: Simulink diagram of the ICE pilot/vehicle/SMC system 
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Figure 12: Bode plots, — , ICE hedge design 
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Figure 17: <t)-tracking task HQ and PIO predictions, ICE failed system 
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Figure 1 8: Outer loop tracking ( a, and p ), ICE, with loop-shaped controller, actuators 

noise, failure at t - 20 sec 
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